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Abstract 

In this paper we investigate the hadron and photon production at forward rapidi- 
ties in pp and pA collisions considering two phenomenological models based on the 
Color Glass Condensate formalism. The predictions of these models for the hadron 
spectra are compared with the dAu and pp RHIC data and the nuclear modifi- 
cation factor RhA is estimated. Our results demonstrate that the analysis of the 
transverse momentum dependence of this factor allows to discriminate between the 
phenomenological models. Furthermore, we predict RhA f° r photon production at 
forward rapidities at RHIC and LHC energies and study the px dependence of the 
photon to pion production ratio. 
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1 Introduction 



Quantum Chromodynamics (QCD) at high energies can be described by an ef- 
fective theory denoted Color Glass Condensate (CGC), which is a many-body 
theory of partons which are weakly coupled albeit non-perturbative due to 
the large number of partons (For reviews see Ref. pQ). Properties of the CGC 
are specified by correlation functions of gluons which evolve with increasing 
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energy. They obey an infinite hierarchy of non-linear evolution equations - the 
Balitsky-JIMWLK hierarchy [2113] . In the mean field approximation, the first 
equation in the hierarchy decouples and it boils down to a single nonlinear 
integro-differential equation: the Balitsky-Kovchegov (BK) equation [21H] . In 
particular, the BK equation determines in the large-iV c limit the evolution 
of the two-point correlation function, which corresponds to scattering ampli- 
tude Af(Y, r) of a dipole off the CGC, with r the dipole size and Y oc In s is 
the rapidity. This quantity encodes information about the hadronic scatter- 
ing and thus about the non-linear and quantum effects in the hadron wave 
function. Although the general solution to the BK equation still is not known, 
approximate solutions have been constructed which separately cover the non- 
linear regime deeply at saturation and the linear regime, where M obeys the 
BFKL or DGLAP equation [5|5f7|IS|9|10] . The transition among these regimes 
is specified by a typical scale, which is energy dependent and is called satu- 
ration scale Q s [Ql oc A a x~ x ]. Basically, the main properties of the solutions 
of the BK equation are: (a) for the interaction of a small dipole (r <C 1/Q S ), 
M{Y = hxl/x, r) « r 2 , implying that this system is weakly interacting; (b) 
for a large dipole (r 3> 1/Q S ), the system is strongly absorbed and therefore 
N(Y,r) « 1. This property is associated to the large density of saturated 
gluons in the hadron wave function. Moreover, the BK equation predicts the 
geometric scaling regime: at small values of x, instead of being a function of a 
priori the two variables r and x, MiY, r) is actually a function of a single vari- 
able r 2 Q 2 s (x) up to inverse dipole sizes significantly larger than the saturation 
scale. This scaling is obvious at r ^> 1/Q S , but it is a non-trivial prediction for 
r < 1/Q S . Furthermore, it breaks down for r <C 1/Q S (leading-twist regime), 
which implies a limited extension for the geometric scaling window. The scal- 
ing behavior is predicted to hold approximately in the range r gs < r < r s , 
where r gs « 1/Q gs (Q gs = Ql/h) and r s « l/Q s - The so-called extended 
scaling region is characterized by the geometric scaling momentum Q gs , which 
grows faster than the saturation scale with x and defines the upper bound in 
transverse momentum of the geometric scaling region. 

The search of signatures for the parton saturation effects has been an ac- 
tive subject of research in the last years. In particular, the geometric scal- 
ing window has been observed in inclusive and diffractive processes at HERA 
|H|ll2|ll3j and the observed [T^ffT5] suppression of high px hadron yields at for- 
ward rapidities in dAu collisions at RHIC had its behavior anticipated on the 
basis of CGC ideas [16J. A current shortcoming of these analyzes comes from 
the non-existence of an exact solution of the non-linear equation in the full 
kinematic range, which implies the construction of phenomenological models 
satisfying the asymptotic behavior which is under theoretical control. Several 
models for the forward dipole cross section have been used in the literature in 
order to fit the HERA and RHIC data [5,6, 7, 8, 9, ID]. In general, the adjoint 
dipole scattering amplitude is modeled in the coordinate space in terms of a 
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simple Glauber-like formula as follows 

N^ipc, r) = 1 — exp 



_l (r 2g2 )7 (^) 



where 7 is the anomalous dimension of the target gluon distribution. The 
fundamental scattering amplitude Mt can also be parameterized as in (CQ), 
with the replacement Q 2 S — > Q 2 s Cf/Ca = 4/9 Q 2 S . Moreover, it is useful to 
assume that the impact parameter dependence of M can be factorized as 
Afp^(x,r,b) = J\f^ t ji(x,r)S(b). The main difference among the distinct phe- 
nomenological models comes from the predicted behavior for the anomalous 
dimension, which determines the transition from the non-linear to the ex- 
tended geometric scaling regimes, as well as from the extended geometric 
scaling to the DGLAP regime. It is the behavior of 7 that determines the fall 
off with increasing pt of the cross section. The current models in the literature 
consider the general form 7 = 7^ + A7, where 7^ is the anomalous dimension at 
the saturation scale and A7 mimics the onset of the geometric scaling region 
and DGLAP regime. One of the basic differences between these models is as- 
sociated to the behavior predicted for A7. While the models proposed in Refs. 
[6117118] assume that A7 depends on terms which violate the geometric scaling, 
i.e. depends separately on r and rapidity Y, the model recently proposed in 
Ref. [10J consider that it is a function of rQ s . In particular, these authors 
demonstrated that the RHIC data for hadron production in dAu collisions for 
all rapidities are compatible with geometric scaling and that geometric scaling 
violations are not observed at RHIC [10]. In contrast, the IIM analysis implies 
that a substantial amount of geometric scaling violations is needed in order to 
accurately describe the ep HERA experimental data [6]. These contrasting re- 
sults motivate the study of other experimental observables directly associated 
to the CGC physics. 

Our goal in this paper is twofold. Firstly, to compare the CGC predictions 
for the ratio RhA [see definition in Eq. ([7]) below] with the RHIC experimen- 
tal data for hadron production. Although the pt dependence of the hadron 
spectra in dAu collisions is quite well described, the overall magnitude of the 
production cross section still is uncertain, which is associated to the fact that 
the CGC calculations have been performed at leading order. Moreover, the 
description of hadron production in pp collisions at forward rapidities using 
CGC physics still is an open question (See e.g. [H]). These two aspects limit 
the predictive power of the CGC formalism for the behavior of the ratio RhA- 
To overcome this, we extend the model proposed in Ref. [TO] for pp colli- 
sions and compare with the recent STAR data for inclusive 7r° production at 
77 = 3.3, 3.8 and 4.0 |15j . It allows to determine the magnitude of the next-to- 
leading order corrections in terms of the .fT-factor necessary to describe the pp 
cross section. A similar calculation is performed for dAu collisions, obtaining 
an identical value for the K-factor at 77 = 4.0. It allows to predict the behavior 
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and normalization of RhA for this rapidity. Assuming that this behavior also 
is present for other rapidities and for charged hadron production, we compare 
the CGC predictions with the BRAHMS data [Hj. Our results demonstrate 
that the study of the pr dependence of the ratio RhA allows to discriminate 
between the distinct phenomenological models. 

Our second goal is to present the predictions of the CGC physics for photon 
production using a model for the scattering dipole amplitude which describes 
quite well the hadron production. As emphasized in Refs. [T81ll9f20j . it is essen- 
tial to consider the electromagnetic probes of the CGC in order to determine 
the dominant physics in the forward region at RHIC and LHC. Distinctly from 
hadron production, there is no hadronization of the final state present in the 
description of the photon production cross section, which implies that it is a 
cleaner probe of the CGC. We estimate the ratio RhA for photon production 
at forward rapidities for RHIC and LHC energies and compare its behavior 
with that predicted for hadrons. Moreover, as a by product, we estimate the 
photon to pion production ratio and study its pr dependence. 

The paper is organized as follows. In the next section (Section [2]) we briefly 
review the hadron and photon production in the Color Glass Condensate for- 
malism and the main characteristics of the distinct parameterizations for the 
dipole amplitude scattering. In Section [3] we define the nuclear modification 
factor RhA and discuss the theoretical expectations for the behavior of this 
ratio. Moreover, we calculate the inclusive 7r° production in pp collisions using 
the CGC formalism and estimate the ratio RhA for pions and charged hadrons. 
Our predictions are compared with the STAR and BRAHMS data. The ratio 
RhA for photon production is estimated in Section H] and the photon to pion 
production ratio is calculated. Finally, in the Section [5] our main conclusions 
are summarized. 



2 Hadron and Photon Production in the Color Glass Condensate 
Formalism 

Lets consider the hadron production at forward rapidity in dAu collisions. 
As pointed in Ref. [2T) . it is a typical example of a dilute-dense process, 
which is an ideal system to study the small- £ components of the Au target 
wave function. In this case the cross section is expressed as a convolution of 
the standard parton distributions for the dilute projectile, the dipole-hadron 
scattering amplitude (which includes the high-density effects) and the par- 
ton fragmentation functions. Basically, the minimum bias invariant yield for 
single-inclusive hadron production in hadron-hadron processes is described in 
the CGC formalism by [51122] 
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Xp 

+ fg/p{xi,PTWA (x 2 , -^-Prj D h/g (~'Pt) 



Xf J \Xf' 

(2) 



where pt, y and xf are the transverse momentum, rapidity and the Feynman- 
x of the produced hadron, respectively. The variable X\ denotes the momen- 
tum fraction of a projectile parton, f(x\,pj) the projectile parton distribu- 
tion functions and D(z,p^) the parton fragmentation functions into hadrons. 
These quantities evolve according to the DGLAP |23j evolution equations and 
respect the momentum sum-rule. In Eq. (j2J), Afp(x,k) and Af^(x,k) are the 
fundamental and adjoint representations of the forward dipole amplitude in 
momentum space, which represent the probability for scattering of a quark 
and a gluon off the nucleus, respectively. Moreover, xp = ^e y and the mo- 
mentum fraction of the target partons is given by X2 = x\e~ 2y (For details see 
e.g. 



The photon production can be evaluated in a similar way pUfTS] , and the 
minimum bias invariant yield can be written in the form [18] 



d 2 NP A^X 



dyd 2 pT (2vr) 2 J xp 



(3) 



where pt and y are now the transverse momentum and rapidity of the pro- 
duced photons. In this equation, D^/ q is the quark-photon fragmentation func- 
tion. Distinctly from hadron production, the rate of photon production only 
depends of the quark content on the projectile hadron. It implies that for the 
region where the gluon contribution can be disregarded in hadron produc- 
tion, the behavior for hadron and production is expected to be similar |18j . 
Furthermore, the two cross sections are dependent on the fundamental dipole 
scattering amplitude, which is a building block of the CGC formalism. There- 
fore, if is constrained for instance in hadron production, the calculation 
of the photon production is straightforward. 

The Eqs. fl2]) and ([3]) are only applicable to forward/backward rapidities in pp 
collisions. On the other hand, in hadron-nucleus collisions at high energies, due 
to the A dependence of the saturation scale, they are expected to also be valid 
for mid-rapidity. It is important to emphasize that the minimum bias cross 
sections discussed in our paper are obtained by impact-parameter averaging 
the inclusive hadron/photon production cross section, which in the CGC for- 
malism depends on the impact parameter only through the saturation scale. 
In Ref. [22] the author discuss two alternatives to implement this calculation, 
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which implies different values for the effective saturation scale (Q 2 S ) in mini- 
mum bias collisions. In what follows we assume that (Q 2 ) = A 1 J^Ql(x /x2) x , 
with A e ff = 18.5(20.0) for dAu (pPb) collisions, in order to compare our 
predictions with those obtained in Refs. [SlfTO] . Moreover, as in Ref. [IT] , we 
assume A e ff = 1 in the pp case. 

The basic input for the calculations of the hadron and photon production 
are the dipole scattering amplitudes Ma an d Mr, which are solutions of the 
Balitsky-JIMWLK hierarchy or the BK evolution equation in mean-field ap- 
proximation. As already explained in the introduction, the general solution to 
the BK equation still is not known, which implies that is necessary to consider 
phenomenological models, based on CGC physics, in order to calculate the ob- 
servables. In what follows we consider two distinct phenomenological models 
constructed to describe the RHIC data: the DHJ model [8] and the recently 
proposed BUW model [10]. In these two models the adjoint dipole scattering 
amplitude in the momentum space is given by 



M a (x,Pt) 



d 2 re ipT - 



oxp i -±(r 2 Q 2 s (x)y^) 



(4) 



where 7 is assumed a function of px rather than r in order to compute the 
Fourier transform more easily. The fundamental scattering amplitude Mp is 
obtained from Ma by the replacement Q 2 — > Q 2 Cf/Ca = 4/9 Q 2 . Moreover, 
in these models it is assumed that j(pt, x) = 7 S + Aj(p T , x), with 7 S = 0.628. 
In the DHJ model, Aj^pt, x) is given by [8] 



Aj DH j(pt,x) = (1 - Is 



\og(p 2 T /Q 2 (x)) 



Xy + d^/y + log{p 2 T /Q 2 {x)) 



(5) 



with y = log(l/x), A = 0.3 and d = 1.2. On the other hand, in the BUW 
model [TH] 



A-f BUW (p T ,x) = (1-7,)— — — — -, (6) 

(oj a — 1) + b 

where u = Pt/Qs{x) and the two free parameters a = 2.82 and b = 168 are 
fitted in order do describe the RHIC data on hadron production. The main 
difference between the parameterizations is the presence of terms in the DHJ 
model which violate the geometric scaling. Distinctly from the BUW model, 
which assumes that A7 satisfies the geometric scaling property, the DHJ one 
predicts that it behaves as \og{p\ / Q 2 (x)) / y for large y and p\ > Q 2 , violating 
the geometric scaling. Another important difference is that the large pt limit 
of 7 — > 1 is approached much faster in the BUW model than in the DHJ one, 
which implies different predictions for the large pt slope of the hadron and 
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photon yield. As shown in [10] , both models describe quite well the dAu RHIC 
data for forward rapidities (y > 2.2), but the DHJ model fails to describe the 
large pt data for smaller rapidities, where the x<i values probed are not very 
small. In next section we extend this analysis for pp collisions and calculate 
the nuclear modification factor RhA for hadron production. 



3 The Nuclear Modification Factor for Hadron Production 

In order to disentangle the nuclear medium effects it is useful to compare 
the data from hadron-nucleus (hA) collisions to proton-proton {pp) using the 
nuclear modification factor RhA defined as: 



where N co u is the number of binary collisions at a given centrality in a hA col- 
lision, which is obtained from a Glauber model calculation (See e.g. Appendix 
I in Ref. |25j). In the absence of nuclear effects, hard processes, as for example 
hadron and photon production at large transverse momentum, scale with the 
number of binary collisions, which implies that RhA = 1 in this case [251126] . 
The rapidity behavior of this ratio is directly related to the description of the 
target. If the target is a dilute system, R h A is expected to grow with rapidity 
because the parton associated to the identified hadron has interacted with a 
greater number of gluons, each contributing with a finite amount of transverse 
momentum. In this case we expect that the ratio assumes a value greater than 
one beyond some value of pt- On the other hand, if the target is a saturated 
system, the ratio is expected to decrease in value for larger rapidities, since the 
CGC physics implies a reduction of the nuclear gluon distribution associated 
to the non-linear evolution. 

Lets present a brief review of the main theoretical expectations for the behav- 
ior of RhA associated to CGC physics (For a detailed review see, e.g, [27J). 
At central rapidities, it is predicted the presence of a Cronin peak, which is 
interpreted as reflecting the classical saturation and is understood as a result 
of Glauber-like multiple scattering off the gluon distribution produced by un- 
correlated valence quarks. The Cronin peak is predicted to disappears after a 
short evolution in rapidity and RhA is suppressed, stabilizing at a small value 
which approaches one asymptotically at large px- The rapid suppression of 
the ratio with increasing rapidity has been interpreted as a consequence of the 
strong difference between the quantum evolution of the nucleus and that of 
the proton. Since for a fixed value of y and pt the proton and nucleus satura- 
tion scales are different, the transverse phase space available for the evolution 
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Fig. 1. Inclusive tt production cross section in pp collisions at RHIC energies. Data 
from STAR collaboration [15]. We assume K{j] = 4.0) = K (rj = 3.8) = 1.4 and 
K(rj = 3.3) = 1.0 for the DHJ (solid line) and BUW (long-dashed line) predictions. 



is larger for the proton than for the nucleus. The amount of suppression is 
estimated as being R^a ~ ^/^loii ■> where 7 is the anomalous dimension which 
depends on the rapidity and transverse momentum. Therefore 7 determines 
the maximal possible suppression of the nuclear modification factor due to 
the saturation effects. As the phenomenological CGC-based models assume 
different behaviors for the anomalous dimension, the analysis of Rha can be 
useful to constrain the QCD dynamics. 

The RHIC data for RdAu [HI confirm the qualitative expectations of CGC 
physics [16J. Although it is a very important evidence for CGC physics, it 
is fundamental to demonstrate the quantitative agreement of the experimen- 
tal data with the CGC predictions. In Ref. [7] the authors have obtained a 
satisfactory description of the BRAHMS data for RdAu assuming that a CGC- 
based description of high-pr hadron production in pp collisions is valid (See 
also [28J ) . This is a strong assumption which should be verified. In principle, 
it is expected that for large rapidities the proton saturation scale assumes a 
large value, which implies a large value for the geometric scaling momentum 
Q gs . Therefore, in this range the extended geometric scaling window becomes 
large and eventually covers the entire regime of particle production, since the 
DGLAP region is cut-off by energy-momentum conservation constraints |17j . 
On the other hand, for mid-rapidity a CGC-based description for pp collisions 
may not be well-justified. Consequently, it is important to test the applicabil- 
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Fig. 2. Inclusive 7r° production cross section in dAu collisions at RHIC energies. 
Data from STAR collaboration [15]. We assume K(r] = 4.0) = 1.4 for the DHJ 
(solid line) and BUW (long-dashed line) predictions. 



ity of the CGC physics in pp collisions at RHIC and verify the rapidity range 
in which this approach can be used. Recently, the STAR collaboration [15] 
has reported the measurements of the production of forward ir° mesons in pp 
and dAu collisions at Wsnn — 200 GeV. These data are ideal to check the 
CGC predictions. In Fig. [T] we compare the DHJ and BUW predictions for 
the minimum bias invariant yield with the STAR pp data. In our calculations 
we use the CTEQ5L parameterization (29] for the parton distribution func- 
tions and the KKP parameterization for the fragmentation functions |30j . As 
in previous calculations [SlfTU] there is one free parameter in our calculation: 
the i^-factor. It is determined in order to obtain the better description of the 
experimental data and is fixed for each rapidity. We can see that the DHJ 
and BUW predictions are almost identical at the two larger values of rapidity, 
but differ in the large pt region for 77 = 3.3, with the DHJ one being larger 
than the data. This trend is similar to observed for dAu collisions in [10J with 
increasing rapidity. 

Some comments are in order here. First, the i^-factors for the different ra- 
pidities were fixed in order to describe the low-p^ data, since they are in 
the extended scaling region where the formalism is expected to be valid (See 
Fig. 2 from Ref. [H]). Second, the .fT-factor necessary to describe the ex- 
perimental data at 77 = 4 is identical that found in Ref. [T7j . where the 
DHJ model was applied to describe pp collision. Finally, we have found that 
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Fig. 3. Nuclear modification ratio RdAu f° r charged hadrons and tt° production at 
RHIC energies. Data from BRAHMS [14] and STAR [15J collaborations. 

K(i] = 3.3) < K(r] = 4.0). This behavior is opposite to the observed when we 
apply the CGC formalism for charged hadron production in dAu collisions. 
At this moment, we were not able to find a reasonable explanation for this 
particular behavior required to describe the pp data at rj = 3.3. 



A similar study can be performed for the 7r° production in dAu collisions. As 
in Refs. PfTU] we assume isospin invariance to obtain the parton distributions 
for a deuteron from those for a proton. In Fig. [2] we present a comparison 
between the DHJ and BUW predictions and the STAR data for 77 = 4. We 
have that the two CGC-based predictions are very similar in this range, as 
already verified in [TU|. An interesting aspect is that the .fT-factor necessary to 
describe the dAu data is identical to that used in the description of the pp data 
at the same rapidity. It implies that the resulting CGC prediction for the ratio 
RdAu at this rapidity would be independent of the ^-factor. Moreover, the pr- 
behavior of this ratio would be a robust prediction of the CGC approaches. 
In Fig. [3] (right panel) we present our predictions for the ratio RdAu for 7r° 
production and i] = 4, where we have assumed that N co u = 7.2 as useful in 
the experimental analysis [H]. We have that the normalization and the pt 
dependence of the experimental data are quite well described by the CGC- 
based predictions. It is a strong evidence for the CGC physics in the forward 
rapidity at RHIC. However, in order to discriminate between the DHJ and 
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BUW predictions we need to consider a larger range of rapidities. 

Motivated by the satisfactory description of the tt° data in pp collisions we 
extend our analysis for charged hadron production at rj = 2.2 and 3.2, where 
we still expect that a CGC calculation is valid. A current shortcoming is that 
there are not experimental data available in literature for charged hadron pro- 
duction at forward rapidities in pp collisions. Therefore, it is not possible to 
constrain the .ff -factor for these cross sections. On the other hand, for dAu 
collisions, the charged hadron spectra were studied by the BRAHMS collab- 
oration [14] . We have calculated the corresponding cross section and verified 
that the BUW model describe quite the data, while the DHJ model fails for 
central rapidity, as already verified in [10]. The basic difference between our 
results and those from [10] is that we have found a X-factor which is two 
times larger than that obtained in [TO], which is directly associated to the 
treatment for the deuteron contribution to the cross section. In our case we 
have assumed that this contribution is normalized by the atomic number. It 
explains the difference by a factor two of our f^-factor at rj = 4 and that 
quoted in [10]. A comment is order here. We have estimated the contributions 
of J\fj?(x, k) and Na(x, k) for the charged hadron cross section in dAu colli- 
sions considering the BUW model and observed that, similarly to the DHJ 
one, the Af^(x, k) contribution determines the large pt behavior of the cross 
section for forward rapidities at RHIC energy, while N"^(x, k) is the relevant 
contribution at mid-rapidity. 

In order to calculate the ratio RdAu for charged hadron production and to com- 
pare with the BRAHMS data [H] we assume that the i^-factor is the same in 
our dAu and pp calculations. This assumption is not trivial: as the saturation 
scale of the nucleus and the proton are distinct, different dynamical effects 
are being probed for a fixed rapidity. Consequently, the normalization of our 
calculations of RdAu for charged hadron can be modified in the future. On the 
other hand, we believe that the px dependence predicted by the CGC physics 
is a robust result which is directly associated to the anomalous dimension 
considered in the distinct phenomenological models. In Fig. [3] (left and mid- 
dle panels) we present our predictions for RdAu in charged hadron production 
using the DHJ and BUW models. We have that the px dependence predicted 
by these models is very distinct. While the DHJ model predict a ratio which 
is basically pt and rj independent, the BUW model predicts a strong pt de- 
pendence, with RdAu increasing almost linearly with pt, approaching one to 
large transverse momentum. Moreover, the BUW model also predicts a ra- 
pidity dependence for the ratio, with the slope increasing at smaller values 
of rapidity. These behaviors are observed in the experimental data. It is im- 
portant to emphasize that both models describe the dAu spectra for r] = 2.2 
and 3.2 as shown in [10] and verified in our calculations. Consequently, the 
Pt dependence of the ratio is directly associated to the distinct predictions 
for the pt spectra in pp collisions, which are different already at rj = 3.3, as 
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Fig. 4. Nuclear modification ratio RdAu for the photon production at RHIC energies. 

verified in Fig. [U The reasonable agreement between the BUW model and 
the experimental data is a strong evidence of the CGC physics. Moreover, it 
indicate that the dipole scattering amplitude satisfies the geometric scaling 
property in the forward RHIC kinematical range. 



4 The Nuclear Modification Factor for Photon Production 



The minimum bias yield for photon production in the CGC formalism can be 
calculated using the Eq. ([3]). The basic input is the fundamental scattering 
amplitude A/jr, which also is present in the calculations of hadron production 
cross sections. In particular, at forward rapidities it determines the behavior of 
this cross section, since the projectile gluon distribution vanishes at x\ — *■ 1. 
In the previous section we have estimated the differential cross section for 
hadron production and obtained a quite well description of the px spectra for 
pp and dAu collision at forward rapidities, which implies that the behavior of 
Mr is reasonably well determined. It allows to obtain reliable predictions for 
the behavior of the photon production cross section. Currently, experimental 
results at RHIC shown that the prompt photon cross section at mid-rapidity 
scale with N co u [31j, which indicate that the nuclear effects are small at rj = 
0. On the other hand, there is not available experimental data for photon 
production at forward rapidities in pp and dAu collisions. We focus our analysis 
in the calculation of the ratio RhA, as defined in the Eq. ([7]) above, at forward 
rapidities. Basically, we calculate the pp and hA minimum bias yields for 
photon production using Eq. ([3]), the CTEQ5L parameterization [29] for the 
parton distribution functions and the GRV parameterization for the quark- 
photon fragmentation function [32]. Similarly to hadron production we assume 
that the K factor is the same for pp and hA collisions. In the particular case of 
dA collisions, we again assume N co u = 7.2 and the isospin symmetry in order 
to calculate the parton distributions of deuteron (For a recent discussion about 
isospin effects in prompt photon production in AA collisions see [33]). 
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Fig. 5. Nuclear modification ratio R p pt for charged hadrons, tt° and photons at LHC 
energies (y/s = 8.8 TeV) and rapidity y=6, considering the BUW model. 



Initially lets calculate the ratio RdAu for photon production at RHIC energies 
(s/snn = 200 GeV) and forward rapidities (y = 2.0, 3.0, 4.0). At smaller 
rapidities, a CGC description for pp collisions is expected to breaks down. In 
Fig. H] we present our predictions for RdAu using the DH J (left panel) and the 
BUW model (right model). There is a large difference between the behaviors 
predicted by the two models. While the DHJ model predicts an almost flat 
ratio, which is pt and y independent, the BUW model predicts that the ratio 
is flat only at very large rapidities, increasing with px at smaller values of 
rapidity. This behavior is similar to that observed for hadron production. 
Consequently, the study of photon production can be an important search 
of information about the behavior of the scattering amplitude and the CGC 
physics. 

A shortcoming for the quantitative understanding of the CGC physics at RHIC 
is associated to the limited phase space in transverse momenta, which implies 
that the transitions expected to occur between the saturation, extended ge- 
ometric scaling and DGLAP regimes are not easily observed. In contrast, at 
LHC energies the available phase space will be much larger even at large ra- 
pidities, allowing to study the different regimes of the QCD at high energies 
in more detail. Here we study the charged hadron, ir° and photon production 
in pp and pPb collisions at ^/s = 8.8 TeV and y = 6 using the BUW model 
and postpone a more detailed analysis for a future publication. For pPb col- 
lisions we assume A e ff = 20 and N co u = 7.4 as quoted in the Table 6 of the 
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Fig. 6. Ratio between the photon and ir° production cross sections at RHIC energies, 
considering two values of rapidity at RHIC (y = 2 and y = 4) and y = 6 at LHC 
energy. 



Appendix I from [25]. In Fig. [5] we present our predictions for the ratio R p pb- 
We can see that the magnitude and dependence is almost identical for the 
different observables. A similar result was obtained in Ref. [31], where it was 
observed almost the same suppression as a function of transverse momentum 
for gluons and heavy quarks. Moreover, we observe that the ratio increases 
with the transverse momentum, as already verified at RHIC. However, the 
ratio is almost one only at px > 10 GeV, which is directly associated to the 
larger window of the extended geometric scaling regime for the proton and 
nucleus at LHC energies. 

Finally, as a by product we calculate the ratio between the photon and hadron 
cross sections. Distinctly from Ref. [35] we focus here in the transverse mo- 
mentum dependence of this ratio at fixed rapidity for dAu and pPb collisions 
at RHIC and LHC energies. It is expected that the i^-factor cancels in this ra- 
tio, which implies that its behavior should not be modified by next-to-leading 
corrections. Following Ref. [18] we focus in the low pt region and forward 
rapidities, where the fragmentation contribution for photon production is ex- 
pected to contributes significantly for the produced photons. In Fig. 6 we 
present our predictions for ratio 7/vr , calculated using the Eqs. (T2]) and ([3]) 
considering the DHJ and BUW models. Due to the distinct phase space avail- 
able for the different rapidities and energies, the curves in the figure finish 
in different points. We can see the DHJ and BUW results are similar, with 
the ratio increasing in the small pt region and saturating at large values of 
the transverse momentum. It means that the ratio 7/71"° is less sensitive to the 
phenomenological model used as input in our calculations and is mainly deter- 
mined by the photon and hadron fragmentation functions. Moreover, the ratio 
increases with the rapidity, as already verified in [35 J, which is associated to 
the fact that the gluon contribution in the projectile hadron diminishes with 
the rapidity. Finally, it is interesting to observe that the distinct predictions 
for the ratio tends to a same value at large pr- 
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5 Conclusions 



The observed suppression of the normalized hadron production in dAu colli- 
sions as compared to pp collisions has been considered an important signature 
of the Color Glass Condensate physics. In the last years several models were 
proposed to describe the hadron spectra in dAu collisions, obtaining a satis- 
factory description of these experimental data. In general, these models have 
been extended for pp collisions in order to calculate the ratio RhA without 
a comparison with the corresponding experimental data. In this paper we 
have, for the first time, estimated the hadron production in pp and dAu col- 
lisions in a same theoretical formalism and compared these predictions with 
the experimental data. The comparison with the STAR data for 7r° produc- 
tion allows to fix the free parameter in our calculations (the .fT-factor) and 
obtain a parameter free prediction for the nuclear modification ratio RhA at 
rj = 4. For other rapidities, there are not available, simultaneously, pp and 
dAu experimental data for hadron production. In order to calculate the ratio 
RhA for these rapidities we have assumed that the A^-factor is the same for 
pp and hA collisions. As discussed before, it is not a trivial assumption. How- 
ever, the predictions for the p^-dependence of the ratio are not affected by 
this choice. The comparison with the experimental data demonstrate that the 
BUW model, which assumes the geometric scaling property, is the adequate 
one for the RHIC kinematical range. 



We also have investigated the photon production at forward rapidities, which 
is considered a cleaner probe of the CGC physics. We demonstrate that the 
behavior of the ratio for photons is similar to verified for hadrons. It implies 
that the study of photon production is a useful search of information about 
the basic building block of the CGC formalism: the fundamental scattering 
amplitude. Finally, as a by product, we have estimated the ratio between the 
photon and hadron cross sections and demonstrated that is not sensitive to 
the phenomenological model used as input in the calculations. 



Acknowledgments 



This work was partially financed by the Brazilian funding agencies CNPq and 
FAPERGS. 



15 



References 



[I] E. Iancu and R. Venugopalan, |arXiv:hep- ph/0303204.; H. Weigert, Prog. Part. 
Nucl. Phys. 55(2005) 461; J. Jalilian-Marian and Y. V. Kovchegov, Prog. Part. 
Nucl. Phys. 56 (2006) 104. 

[2] I. I. Balitsky, Nucl. Phys. B463 (1996) 99, Phys. Rev. Lett. 81 (1998) 2024, 
Phys. Rev. D 60 (1999) 014020, Phys. Lett. B 518 (2001) 235; I.I. Balitsky and 
A.V. Belitsky, Nucl. Phys. B629 (2002) 290. 

[3] E. Iancu, A. Leonidov and L. McLerran, Nucl.Phys. A692 (2001) 583; E. 
Ferreiro, E. Iancu, A. Leonidov and L. McLerran, Nucl. Phys. A701 (2002) 
489; J. Jalilian-Marian, A. Kovner, L. McLerran and H. Weigert, Phys. Rev. D 
55 (1997) 5414; J. Jalilian-Marian, A. Kovner and H. Weigert, Phys. Rev. D 
59 (1999) 014014, ibid. 59 (1999) 014015, ibid. 59 (1999) 034007 ; A. Kovner, 
J. Guilherme Milhano and H. Weigert, Phys. Rev. D 62 (2000) 114005; H. 
Weigert, Nucl. Phys. A703 (2002) 823. 

[4] Y.V. Kovchegov, Phys. Rev. D 60 (1999) 034008, ibid. 61 (2000) 074018. 

[5] K. Golec-Biernat and M. Wusthoff, Phys. Rev. D 59 (1999) 014017, ibid. D60 
(1999) 114023; J. Bartels, K. Golec-Biernat, H. Kowalski, Phys. Rev. D 66 
(2002) 014001; H. Kowalski and D. Teaney, Phys. Rev. D 68 (2003) 114005; 
H. Kowalski, L. Motyka and G. Watt, Phys. Rev. D 74(2006) 074016; K. Golec- 
Biernat and S. Sapeta, Phys. Rev. D 74 (2006) 054032; J. T. de Santana Amaral, 
M. B. Gay Ducati, M. A. Betemps and G. Soyez, Phys. Rev. D 76 (2007) 
094018; C. Marquet, R. Peschanski and G. Soyez, Phys. Rev. D 76 (2007) 
034011; G. Soyez, Phys. Lett. B 655 (2007) 32. 

[6] E. Iancu, K. Itakura, S. Munier, Phys. Lett. B590 (2004) 199. 

[7] D. Kharzeev, Y.V. Kovchegov and K. Tuchin, Phys. Lett. B599 (2004) 23. 

[8] A. Dumitru, A. Hayashigaki and J. Jalilian-Marian, Nucl. Phys. A 765 (2006) 
464; Nucl. Phys. A 770 (2006) 57. 

[9] V. P. Goncalves, M. S. Kugeratski, M. V. T. Machado and F. S. Navarra, Phys. 
Lett. B 643 (2006) 273. 

[10] D. Boer, A. Utermann, E. Wessels, Phys. Rev. D 77 (2008) 054014 

[II] A. M. Stasto, K. Golec-Biernat and J. Kwiecihski, Phys. Rev. Lett. 86 (2001) 
596. 

[12] C. Marquet and L. Schoeffel, Phys. Lett. B639 (2006) 471. 

[13] V. P. Goncalves and M. V. T. Machado, Phys. Rev. Lett. 91 (2003) 202002; 
JHEP 0704(2007) 028. 

[14] I. Arsene et al. [BRAHMS Collaboration], Phys. Rev. Lett. 91 (2003) 072305; 
Phys. Rev. Lett. 93 (2004) 242303; Phys. Rev. Lett. 94 (2005) 032301. 



16 



[15] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 97 (2006) 152302 

[16] R. Baier, A. Kovner and U. A. Wiedemann, Phys. Rev. D 68 (2003) 
054009; J. Jalilian-Marian, Y. Nara and R. Venugopalan, Phys. Lett. B 577 
(2003) 54; D. Kharzeev, Y. V. Kovchegov and K. Tuchin, Phys. Rev. D 68 
(2003) 094013; J. L. Albacete, N. Armesto, A. Kovner, C. A. Salgado and 
U. A. Wiedemann, Phys. Rev. Lett. 92 (2004) 082001; E. Iancu, K. Itakura 
and D. N. Triantafyllopoulos, Nucl. Phys. A 742 (2004) 182. 

[17] D. Boer, A. Dumitru and A. Hayashigaki, Phys. Rev. D 74 (2006) 074018 

[18] J. Jalilian-Marian, Nucl. Phys. A 753 (2005) 307 

[19] M. A. Betemps and M. B. Gay Ducati, Phys. Rev. D 70 (2004) 116005; Phys. 
Lett. B 636 (2006) 46. 

[20] F. Gelis and J. Jalilian-Marian, Phys. Rev. D 76 (2007) 074015 

[21] E. Iancu, C. Marquet and G. Soyez, Nucl. Phys. A 780 (2006) 52 

[22] A. Hayashigaki, Nucl. Phys. A 775 (2006) 51 

[23] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438; G. Altarelli 
and G. Parisi, Nucl. Phys. B126 (1977) 298; Yu.L. Dokshitzer, Sov. Phys. JETP 
46 (1977) 641. 

[24] F. Gelis and J. Jalilian-Marian, Phys. Rev. D 66 (2002) 014021 

[25] F. Arleo et al, "Photon physics in heavy ion collisions at the LHC," 
|arXiv:hep-ph/031113T| 

[26] A. Accardi et al, "Hard probes in heavy ion collisions at the LHC: PDFs, 
shadowing and p A collisions," |arXiv:hep-ph /0308248. 

[27] E. Iancu, AIP Conf. Proc. 739 (2005) 362 [arXiv:hep-ph/0408228] . 

[28] N. Armesto, C. A. Salgado and U. A. Wiedemann, Phys. Rev. Lett. 94 (2005) 
022002 

[29] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12 (2000) 375. 
[30] B. A. Kniehl, G. Kramer and B. Potter, Nucl. Phys. B 582 (2000) 514. 
[31] P. Stankus, Ann. Rev. Nucl. Part. Sci. 55 (2005) 517. 

[32] M. Gluck, E. Reya and A. Vogt, Phys. Rev. D 48 (1993) 116 [Erratum-ibid. D 
51 (1995) 1427]. 

[33] S. Turbide, C. Gale, E. Frodermann and U. Heinz, Phys. Rev. C 77 (2008) 
024909. 

[34] K. Tuchin, Nucl. Phys. A 798 (2008) 61. 

[35] J. Jalilian-Marian, Nucl. Phys. A 806 (2008) 305. 



17 



